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1.0 INTRODUCTION 


Many studies, such as that by Chun, et.al, (Ref. 1), have been n'.ade to 
analyze jet engine noise. In general, these studies are conducted in a static 
mode, and predictions are made to reflect an actual flight situation. The re- 
sult is usually expressed as a value of EPNdB (Effective Perceived Noise in 
decibels), where much of the detailed information is "lost" due to the re- 
quirements for computing EPNdB. 

In the analysis package which follows, the intention is to convert flight 
data to equivalent static data using current prediction methods. Tones are 
clearly distinguishable from broadband noise since narrow-band analysis is em- 
ployed. Eventually, such narrow-band analysis of flight data is expected to 
result in criteria for ground tests, which are easier and less costly to per- 
form than flight tests. 

The package presented is one which encompasses several problems associa- 
ted with acoustical analysis of a moving source with respect to a stationary 
observer. The nonstationarity of the data causes difficulty in applying con- 
ventional time series analysis. Propagation effects influence all recorded 
data and must be accounted for. Also, the short integration time for each re- 
cording microphone requires some type of signal enhancement to increase 
accuracy of the data levels. 


2.0 EXPERIMENTAL DESIGN 


The development of this data reduction package was initiated to study jet 
engine flight noise where fan tones were radiating from the engine inlet (see 
Fig. 1). The aircraft's flight track was such that the noise source was flown 
directly over an array of microphones at constant altitude and velocity. 


1 





Figure 1. - General test design 



The aircraft's position, in general, can vary from the intended flight path. 
Those variations were recorded by the use of a laser radar. The coordinates 
of the source (X', Y' , 7') as shown in Figure 2, are defined by the displace- 
ments from an axis system (X, Y, and Z) whose origin is the first microphone 
In the array* These coordinates were obtained from a spherical to Cartesian 
transformation of the radar data followed by a translation of the origin to 
trie first microphone. 

Weather information was obtained with a specially implemented balloon 
stationed near the microphone array. Temperature, barometric pressure, rela- 
tive humidity and wind speed were recorded at various altitudes from ground 
level to the altitude of the aircraft. All but wind speed are used in the 
propagation corrections. 


3.0 STATIONARITY VERSUS NONSTATIONARITY 


Many studies have been completed on the farfield noise levels of a sta- 
tionary jet engine (i.e., via static testing) and the methodology for inter- 
preting the noise characteristics are well understood. When the noise source 
is moving, however, the farfield noise is not as predictable as in the static 
mode since the forward speed effects (the motion effects) on the noise are not 
well understood (Ref. 1). In analyzing aircraft flyover data, it is desirable 
to obtain accurate spectra at various aircraft positions to determine noise 
levels generated by the source and to offer some understanding as to the for- 
ward speed effects. 

A main area of concern in analysis of flyover data is its nonstation- 
arity. Data whose statistical properties vary with the passage of time (Ref. 
2) are known as nonstationary. Frequency analysis, or time series analysis, 
was developed to handle random, stationary data. Hence, some difficulty is 
encountered in determining the frequency content of nonstationary data utiliz- 
ing conventional data reduction techniques. 
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3^1 Iitipl i cation of the Doppler Effect 


■»» 


The Doppler effect is of major concern as the apparent frequency changes 
throughout the movement of the source past the observer. For any microphone, 
the recorded frequency is equal to the actual generated frequency only when 
the aircraft is positioned directly over the microphone (0*90®), Thus, a 
method of analysis must be employed which assumes that a flyover noise data 
set meets the criteria for use of conventional time series analysis whose 
algorithms are generally based upon use of the Fast Fourier Transfonn (FFT). 

The problem of nonstationarity in flyover analysis has been studied by 
such people as E.P. McDaid and L. Maestrello (Ref. 3) who found that if the 
directionality of the source is not tai^t^n into consideration, the effect of 
nonstationarity is negligible for most practical cases. In general, if the 
aircraft is considered to move discretely along its path from point to point, 
one may accept a small increment of time during which the data is relatively 
or locally stationary. According to J. S, Dendat and A, G. Piersol (Ref. 2) 
this assumption is acceptable providing the statistical properties within this 
increment do not change and hence time series analysis may be employed. 
Specifically, as has been shown with vibration data, e.g., vibration of a 
spacecraft during launch (Ref. 2), data may be considered locally stationary 
over those small increments if the data has normal and Chi -squared 
distributions. 


3.2 Ensemble Averaging 


Choosing a small increment of time over which to apply time series analy 
sis inhibits the ability to ensemble average to obtain frequency content in- 
formation. Yet, for most algorithms, it is desirable to average to obtain 
reasonable statistical accuracy. With a linear array of microphones placed 
along the flight path, each, theoretically, will record identical spectra at 


- * a Wli o 


OF 


different recording times, barring any transient occurrences with the source 
production or any atmospheric disturbances. Referring tn Moure 3, that time 
«;enaration At^* is 


At. « d./V 

where d^. is the distance between each microphone and the refe.rence micro- 
phone, and 1 is the veircity d-5»'ect1y the m1c«'ontione array. Spectra at 

corresponding angles, n, mau he av/eraoed to I'^ad fn a roBulting ensemble 
averaged spectra. Thus, the total nunbon nf a'^nranes I. for •^he resultant 
power spectra is the number M of fb'T's averargri m th^ t-Sw series analysis 
times the number M of microphones averaned, o*’ The final power 

spectral density RSllj^ may be ejiprossed as 

N 

PSD. « V psn. 
i-4 ^ 

where PSD^. is the power spectra for micf'ophnne i calculated over M 
ensemble averages. 


^ Pet e rm i nat i on of t h e f't jtnjy^ i^of Averages 

As was previously stated, one would like to acgui'^o some specified level 
of accuracy for each spectral estimate. In general, the resulting power 
spectral density accuracy increases as the number o^ degrees of freedom MD 
increases. For PSD calcuation wia the direct method (Ref. 4), 


NP ?,l 
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Figure 3, - Microphone position and < 
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There are many ways to numerically define the error associated with power 
spectral density estimates generated by conventional time series analysts 
algorithms. (All those described below apply to PSn's calculated by the 
direct method (Ref. 4), not the Blackman-Tukey method.) The first expression 
for the error e of estimation is 


® /T / TT 

This states that the error e decreases as NO Increases. Once NO reaches 
100, or the number of ensemble averages Is 50, little accuracy Is gained, 

A second method results in defining confidence Intervals based upon a 
Chi-squared distribution (Ref. 4). A percentage of confidence, or the percent 
probability that the measured mean square pressure spectrum accurately 
represents the true mean square pressuie spectrum, may be chosen. Depending 
on the number of ensemble averages, a confidence interval at that probability 
is calci'lated* For example, at a 90% confidence level and 5 ensemble averages 
(10 degrees of freedom), the confidence interval is from -4.1 dB to I?.. 9 dB, 

To achieve accuracy within ±1 dB at 90% confidence, 40 averages would be 
needed (actual confidence interval 1s -1.2 dB to 1.0 dB), 

Still another method was studied by K. Rao and J. Preisser (Ref. 6). The 
estimated and asymptotic variances were compared to determine the number of 
averages necessary to produce an adequate spectra. To achieve a reasonable 
normalized random error in percent, 


L w (Bp)^W 
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where W Is the weighting factor for the data window applied to the time do- 
main data when calculating the Discrete Fourl r Transform, Is a measure 
of the convergence of mean spectra. If Is 10, which corresponds to a 9055 
confidence, and the Hann data window .s applied (W ^ 3/B), 

I. w (10^*) t 3/8 « 40 averages 

From the error measurements presented above, a total of 40 averages 
appears to be sufficient to result In a satisfactory spectral representation 
for most flyover noise data. One must remember, however, that the data must 
be relatively stationary over the time interval corresponding to M averages 
or blocks. A block is the time segment over which the Fourier Transform is 
applied. Hence, a case where 5 Fourier Transforms (M«5) per microphone and B 
microphones (N«8) are averaged could be utilized if the stated criteria are 
met . 


4.0 PROPAGATION EFFECTS AND RACKGROUND 
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To obtain an accurate spectral representation of noise data, it Is Impor- 
tant to account for all physical phenomenon present. In this section, a brief 
overview of the propagation effects and background is given. The methods em- 
ployed in the flyover analysis package are discussed and the equations for 
their calculation are given. 


^ 1 The Doppler Effect 


The Doppler effect, or the apparent change in frequency due to the rela- 
tive motion of the source to the observer, is perhaps the most well known of 
the propagation effects. In the case of narrow-band flyover analysis, it can- 
not be ignored. 


. * « 


i. 


As the aircraft passes over a microphone at some average velocity Y, the 
observed frequency f^ is related to the actual source frequency f^ by 

fs “ fo 

where M„ is the Mach number (M. *= Y/c where c is the speed of sound in 
c c 

the medium through which the wave travels) and e is the angle previously 
defined. 


^•2 Convective Amplification 


It has been well established that an acoustic signal is amplified due to 
the motion of the source. Various mathematical expressions for 1 effect 
exist depending on the type of source in motion. 

A simple model which has been used up until very recently incorporates a 
small pulsating sphere represented by a convective monopole. A. Dowling 
(Ref. 7) states that this model is not accurate as motion introduces addi- 
tional coupled monopolas whose effects lead to convective features previously 
not predicted. Convective amplification does depend upon the geometry of the 
source as v/ell. 


At this point, no representation which encompasses all that is discussed 
by Dowling exists. A frequently used expression relating the source pressure 
Pg to the observed pressure Pq is 


(1 - M^cosO) 


2n+2 
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where n Indicates the type of noise source. 

0 Monopole 

n « 1 Dipole 

P, Quadrapole 

The expression (1 - M^coso) is the same as was defined in the Doppler effect. 
The difference in sound pressure level between static and flight cases aSPL^ 
is 


ASPLj = 20{2n+2) lopio (1 - M^cose) 

V 


4.3 Inverse Square Law 


The inverse square law describes the effect of the intensity of a signal 
falling off as 1/r^ where r is the radial distance from the source to the 
observer. In other words, if n and rz correspond to two points on a ray 
emanating from a source, the respective acoustic pressures are related by 


Pz^ = 


P]/- 


(£ 1)2 


The difference in sound pressure level ASPL 2 may then be expressed as 

4SPI.J = -20 (log 10 (f^)) 
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Atmospheric Absorpti on 


Sound absorption in still air leads to an attenuation of the wave as it 
passes through the atmospheric medium. Atmospheric absorption has been 
studied quite extensively. For example, C. M. Harris (Ref. 8) defined the co- 
efficient of absorption a under controlled conditions for various values of 
relative humidity, temperature, and frequency; M. Greenspan (Ref, 9) studied 
the rotational relaxation of nitrogen, oxygen, and air; K. S. Chun, et.al. 
(Ref. 1) offer a simplified calculation of a. All are attempts to quantize 
the total atmospheric absorption into thermal and viscous effects (called 
classical absorption) and rotational and vibrational relaxation effects. Vi- 
brational relaxation is primarily due to both nitrogen and oxygen relaxation. 

The absorption coefficient is a composite of classical absorption 
rotational relaxation “pot* vibrational relaxation of nitrogen and 
oxygen, “vib,0 ''espectively , or 

“ = “CL-" “rot+ “vib.N-" “vib,0 


F. D. Shields and H, F. Bass (Ref. 10) have combined these coefficients to 
provide a thorough method of calculating the absorption coefficient in terms 
of dB/meter which can easily be applied to sound pressure level data. The 
following is an outline of their development. 

Given the barometric pressure P, temperature T, and relative humidity 
RH, at any frequency f, the following procedure may be employed. 
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Calculate the partial pressure of saturated water vapor In N/m^ by 

1»9lo(Psat'''’o) ' 10-79586 [1 - (T^j^/T)] 

-5.02808 logio(T/Toi) + 1.50474 x io~'*(l-10-8*29692C(T/Toi)-l]j 
+0*42873 X 10-3{10‘* 76955Cl-(Toi/T)]^j^ -2.2195983 

where = reference pressure of 1.013 x 10^ N 

Toi = PJS.ie^K 

Calculate the absolute humidity H in % 

M = "H C’sat'''’o>/''/'’o> 

Calculate the relaxation frequency of oxygen and nitrogen, f^ q 
and by 

f n = (P/Pfh [24 + 4.41 X 10+^* H[(0.05 + H)/(0.391 + H)]] 

P jU 0 

^r.N = '”‘9 (-5-l«[{VTg)-l^^-l])] 

where = reference temperature of 293.15°K. 

Calculate the absorption coefficient a(f) in dB/m 
a(f) = 8.686(T/T^)'^^-i:fV{P/Pn)] x 

(1.84 X 10-" + 2.18 X 10-‘'(T/T^)-1(P/P^)(2239/T)2x 
[ exp (-2239/T)]/[f^^n + (fVf^^g)] 


« « • 




+ 8.16 X 10-‘'(T/T|j)-'(P/Pj|)( 3352/T)2 X 
[exp (-3352/T)]/[f^^„ + (fVfp^N)]) 

Once the atmospheric absorption coefficient is calculated, it may be applied 
to atmospheric layers, each with corresponding P, T, and RH value (see Fig. 
4). Over N layers {n) the attenuation of the sound wave in terms of dB, 
ASPL 3 , is 

N 

ASPL 3 = I a.(f) . r. 
i=i 

at each frequency f, where r.j = h/sine. 


4.5 Ground Impedance 

The levels recorded by a microphone above the ground include energy which 
has been reflected by the surface. This additional ini isity must be sub- 
tracted from the observed sound pressure level values to obtain a free field 
value. 

Consider Figure 5 in which both a reflected wave and a direct wave from 
the source are recorded by the microphone. As explained by Pao et.al. (Ref. 
11 ), the ground factor or ratio of the free field mean square pressure to the 
mean square pressure with ground effect, when the surface is considered to be 
acoustically hard is 

G = 2 + cos(kAr) 

where k is the wave number , a = 0.01, and Ar is the difference 
between the reflected path and the direct path. When the noise is averaged 
over finite frequency bandwidths, 
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Microphone 




absorption layered model 
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if h «Z and h«RE 
Figure 5. - Ground impedance model 
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B = 2 + 2 cos(k^r) 


2tt Af ^ 

where Ak » ~ • -g^ (where Af is the bandwidth of interest), and is “ 
(where f^ is the central frequency about the band.)* The resulting 
difference in sound pressure level ASPL^ is 


ASPL 4 == -10 log (2 + 2 ' cos (k,.r) ) 


4.6 Background Subtraction 


It may be desirable to subtract background noise from a calculated sound 
pressure level. This is easily accomplished by comparing a background 
spectrum at some value of 0 to a data spectrum at the same angle. The 
spectral values must be compared at each value of frequency. 


In general, to subtract sound pressure level values at some frequency f|^, 

( S^LfkB"SP*-fkm ) 

SPLf^C = + 10 loqio (1 - 10 ) 


where SPL^I^^ is the background subtracted resultant SPL, SPL^I^^ is the 
measured or data SPL value, and SPL^j^g is the background SPL value. 


For a signal to noise ratio e(f|^), where e(f|^) = SPL^I^^-SPL^I^g , greater 
than 10 dB, a negligible correction is required as the background level is 
considerably smaller than the level of interest. A signal to noise ratio less 
than or equal to 3 dB implies that the background and data Tevels are very 
close, i.e., background only. If the data consistently shows a small signal 
to noise ratio throughout all frequencies, some question may be raised as to 
the validity of either the background spectra or the data spectra. 
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5.0 FLYOVER ANALYSIS DATA REDUCTION PROCEDURES 


As was mentioned earlier, for flyo'.^r analysis narrow-band spectra and 
directivities are the desired output. The directivity is often used as a 
characteristic measurement of a jet engine, and sound pressure levels are used 
for more extensive analysis. The procedures employed to yield these results 
are six-fold. 


5 ^ Analoq-to-Digital to Engineering Units Tape 


The first procedure results in an engineering units tape, for example, a 
tape whose data channel units are N/m^, which places the information in a 
readable form for the performance of all succeeding functions (see Fig. 6). 

The first step involves analog-to-digital eonverslon by available transcrip- 
tion methods. Care must be taken in this step. Before digitizing, considera- 
tion must be given to the maximum frequency of interest and the frequency re- 
solution (bandwidth) to be desired for spectral analysis. Once the maximum 
frequency of interest is known, it is standard practice to low-pass filter at 
that frequency to avoid aliasing, or folding, in the time series analysis re- 
sults. For most hardware systems, a rate digitization, or sample rate SR, 
is required to be 2.5 times the selected cutoff frequency or greater to avoid 
biasing when filtering the data. All time series analyses can yield results 
up to 1/2 • SR which is known as the Nyquist frequency. 


In conjunction with the selection of the maximum frequency of interest 
and hence the sample rate, a frequency resolution must be chosen. The band- 
width BW is 

BW 

uw 2.NPTS 

where NPTS is the number of output points from a time series analysis pro- 
gram. In the direct method of power spectral density computation 
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Analog tape — ^ Digital tape 

(l) Engineering units tape 
(D OVIBSH 

I 

Time “Shifted tape 
Radar file ► AAP (D 

I 

Raw SPL tape 

0 MANDATA 

Averaged raw SPL 

Weather FLYOVER 

Corrected -average SPL 
(' 

© SPLTHTC 
Directivity 

Figure 6. - Data reduction flowchart 
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NPTS is one-half the nuniber of points over which the Fourier Transform is 
applied, i.e., 1/2 the nuinbe*' of points chosen per block. One must reinember 
that over the time interval correspondinq to M blocks, the data should be 
relatively stationary. In other words, the time increment of assumed local 
stationarity tp is 

tp « (2*M*NPTS)/SR 

Hence, maximum frequency, bandwidth, and time of local stationarity must all 
be examined prior to the selection of SR. 

The second step in generating the engineering units tape is to apply the 
proper gains and sensitivities to each recorded data channel. The sensitivi- 
ties are found by recording and digitizing known calibration signals through 
each data channel which yields a linear relationship between counts and engi- 
neering units. Gains are tabulated for each microphone and for each flyover 
which is made. Once these are appl to the digitized data, flyover noise 
analysis may begin. 


5 . t Time Shifting 


To be able to average microphone sound pressure levels, each microphone 
must be shifted in time to appear to be located at the same reference posi- 
tion. This is accomplished by the program called OVIBSH {Appendix A) by 
matrix manipulation (Step 2 of Figure 6), The velocity 7 is extracted from 
radar information and used to calculate the number of points each microphone 
is to be shifted, NSHIFT| , where 

SR*d,. 

NSHIFT. = 

’ ? 

and d^- is the distance from microphone i to the reference microphone or 
reference position. 
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5 ^ ^« aw Sound Prfisy.u r M.^oT^ 


It is at this point, Stop ? of Fiquro 6, that the determination of the 
noise levels begins. A tiinp serins analysis proqram called the Acoustics 
Analysis Program (Ref. 0 and Appendix R) is iitilixed to determine the raw 
sound pressure levels, or the souni pressure levels of tirne«shifted engineer- 
ing units data for each microphone and each selected value of fi. Averaged 
and corrected SPL's and the directivities are calr.ulated from the raw SPU's. 

The Acoustics Analysis Program employs the direct method of computation 
for the power spectral densities, however, that is not a requirement of the 
flyover analysis package. Many time series analysis programs exist which uti- 
lize the Blackman-Tukey method, i.e., the power spectral density is calculated 
from the data set's autocorrelation. Some niffere*n.*'s between the methods 
(Ref. I?.) should be considered to generate cnuparable output, 

^ ^ Averaged gP L 


To obtain corrected sound pressure levels, it is first necessary to 
average the raw sound pressure levels of the goometrically similar microphones 
at each selected angle as discussed in Section 3, This average ms-ist be accom- 
plished by dealing with units of (pressure)'', P''*, or power. More specif- 
ically, the average pov/er P^(f) at frequency f is 

SPlj(f ) 


where N is the number of microphones to he averaged and REF is a reference 

P O 

for dB conversion, (REF = ?. x 10”'’N/m*). This average power may be converted 
back to an average SPL value, SPLg(f) by 
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SPLg(f) ^ JO loqio( 



) 


This Is done by pro^^ram MANOATA (Apppndix C) which ir. Step 4 nf Ftqure 5, 


5 S Corrected Average SPL 


Step 5 of Figure 6 Is the application of program FLYOVER (Appendix D) to 
the averaged sound pressure levels. It corrects tho spectra for Instrumen- 
tation effects and propagation effects. The output is then a realistic pic- 
ture of the source generated noise levels, if background is considered to be 
negligible, i.e., e(f|^) 10 dB for all values nf 

The first correction to be applied is that of instrumentation. It is 
composed of 1) pressure response, ?.) diffraction, and 3) windscreen correc- 
tions. In general, these corrections are frequency dependent and are to be 
added to the observed sound pressure level. They are functions of the type of 
microphone and the angles of acoustic incidence. 

The propagation effects are applied in the following order: 

1) Convective Amplification 

2) Inverse Square Law 

3) Atmospheric Absorption 

4) Ground Impedance 

5) Doppler Frequency Shift 

The order of application is not important with the exception of the Doppler 
effect. Some of the earlier corrections, i.e., atmospheric absorption and 
ground impedance, are frequency donendont and iitili:jo the observed frequency 
for calculation. 


5, <5 Directivity 


The final procedure In Roure 6, Stop i5 involves the calculation of the 
directivity. Program SPLTHTC (Appendix E) determines peak value at a selected 
frequency band fiAf where f is the frequency of interest and Af Is a 
factor which allows for small variations in f from one angle to another as 
Implementation of the Doppler shift does not yield identical values of fre- 
quency for each o. 

Program SPLTHTC also has the capability to subtract the background direc- 
■(•ivlty bv the method discussed In Section 4,6. This results in a true repre- 
sentation of the source's directivity. 

The final option to SPLTHTC is to sum the two largest values within the 
band f±Af. This is necessary to account for a spreading of the peak value to 
two frequency values. The phenomenon is caused by reflection and the fact 
that the Fourier transform is applied to a discrete interval. 

6.0 APPLICATION 


Up to this point, no sample data have been presented. It is the Intent 
of this section to aid in understanding voe data reduction techniques by pre- 
senting an example of flyover analysis (Ref. 13). 

6.1 Design 


For this test, a monotone source of 4000 Hertz was mounted on the wing of 
an aircraft. The source was flown over a microphone array co’Asistirig of 10 
microphones located 30 feet apart and placed 30 feet above the runway (see 
Figr 7). Atmospheric data were recorded by a weather balloon during each 
flight. The aircraft's position was recorded by radar and, in general, was 
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Figure 7. - Microphone array for experifimntal design 
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300 feet above the runway at a velocity of 200 ft/second. The analog tape 
which recorded the pressures of all ten microphones was digitized at 50000 
samples per second. This sample rate was chosen to meet the criteria presen- 
ted in Section 5.1. The maximum frequency of interest was 20000 Hertz and a 
bandwidth of 100 Hertz or less was requested. To be able to average 8 micro- 
phones and 5 transforms for a random error of 10 (see Section 3.3), a block 
size (the number of points over which to compute the Fast Fourier transform) 
of 512 points was chosen. This results in 256 output points, and, 

BW 1^^ = 97.656 Hertz 


and 

^ _ 2*5*256 _ n 0519 cpf. 

^r nrfKTOTJ” sec. 

which corresponds to an aircraft displacement of 10.24 feet along its flight 
path. This is a relatively small distance compared to the total recorded X 
displacement which is approximately 4000 ft. Therefore, it is considered a 
discrete increment for time series analysis purposes. 

6.2 Results 


Given the information above, the microphones must be shifted by a number 
of points equal to 


NSHIFT^. 


(i-l){30 ft) *50000/ sec 
2W ff /f?c 


7500 (i-1) pts. 


where i is the microphone number and microphone 1 is the reference micro- 
phone. (The shifting procedure results in 67,500 fewer digital points and 
should be taken into consideration when determining the time interval for 
di gitization. ) . 
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At this point, raw sound pressure levels are generated by the Acoustics 
Analysis Program. Figures 8 (A through H) show raw sound pressure levels of 
microphones 1 through 8 at an angle of 80". Note that radar information has 
been incorporated and is displayed on the plots. Raw sound pressure levels of 
this type were generated for all microphones over 5 blocks of 512 points at 
angles of 20® through 110° at 5® increments. 

The next step is to average the microphones at each selected angle. 

Figure 9 is the averaged sound pressure level resulting from the spectra seen 
in Figure 8. The averaging process is done in terms of pressure squared. 

Note the averaged spectra's relative smoothness when compared to the spectra 
seen in Figure 8, which is due to the employment of ensemble averaging. 

Program FLYOVER (Step 5 of Figure 6) is now applied to each averaged 
sound pressure level. Input required includes the radar information present 
on each average spectra file, and instrumentation corrections, and weather 
data. Samples of the latter two can be seen in Tables 1 and 2. The averaged 
and corrected sound pressure level results and an example can be seen in 
Figure 10. This is the same data shown in Figures 8 and 9. Note that the 
strongest signal occurs at 4000 Hertz which is the frequency generated by the 
source. 

Once the averaged and corrected sound pressure levels for all selected 
angles have been computed, the directivity may be determined. Choosing 4000 
Hertz as a frequency of interest. Figure 11 results, A value of 200 Hertz 
was given as the band over y/hich to determine the sum of the two highest 
values in program SPLTHTC. The summing is done as the peak is spread over two 
frequency values (see Figure 10). Note that the directivity of the monotone 
source is a constant which is the expected result. 
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(a) Microphone 1 


Figure 8« " Raw sound pressure levels 
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(b) Microphone 2 


Figure 8, - Continued 
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(c) Microphone 3 


Figure 8. - Continued 
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(d) Microphone 4 


Figure 8. - Continued 
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(e) Microphone 5 


Figure 8 
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(f) Microphone 6 


Figure 8. - Continued 
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(g) Microphone 7 


Figure 8, ~ Continued 
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(h) Microphone 8 


Figure 8. - Continued 
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Figure 9. - Averaged sound pressure level 
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Figure 10. - Averaged and corrected sound pressure level 
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TABLE I 


EXAMPLE OF INSTRUMENTATION CORRECTIONS ADDED TO 
MEASURED SOUND PRESSURE LEVELS 


FRQUtNCY 

PRESSURE RESPONSE 

DIFFRACTION 

UINUSCREEN 

(kHz) 

(dB) 

(it) 

(dB) 

1.0 

0.0 

+ .1 

0. 

1.1 

0.0 

+.1 

0. 

1.2 

0.0 

+ . 1 

0. 

1.3 

0.0 

+ .1 

0. 

1.4 

0.0 

+ .1 

0. 

1.5 

0.0 

+ .1 

0. 

1.6 

0.0 

+ .1 

0.^ 

1.7 

0.0 

+ .1 

o' 

1.8 

0.0 

+ .1 

0. 

1.9 

0.0 

+ .1 

0. 

2.0 

o.n 

rt 1 

0. 

2.1 

0.0 

■^ m 2 

-0.3 

2.2 

0.0 

+ .2 

-0.3 

2.3 

0..0 

+ .2 

-0.3 

2.4 

n.o 

+ . 2 

0.0 

2.5 

0.0 

+ .2 

0.0 

2.6 

0.0 

•^.2 

-0.5 

2.7 

0.0 

+ .2 

-0.5 

2.8 

0.0 

+ .2 

-0.5 

2.9 

0.0 

•*-.2 

-0.5 

3.0 

0,0 

+ .2 

-0.5 
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TABLf: 2 


EXAMPLE OF WEATHER DATA NECESSARY FOR PROPAGATION CORRECTION INPUTS 


ALTIiUDE 

TEMPERATURE 

RELATIVE HUMIDITY 

PRESSURE 

■■ 

.{“Cl. 

(%} 

(mbars) 

0 

12.4 

50.0 

1013.0 

10 

12.2 

/in.o 

1011.9 

20 

11. B 

48. B 

1010.7 

30 

11.5 

49.0 

1009.6 

40 

11.0 

51.1 

1008.4 

50 

10.6 

52,0 

1007.4 

60 

10.6 

53,0 

1006.2 

70 

10.6 

53.0 

1005.1 

80 

10.7 

53.0 

1003.9 

90 

10.7 

54.0 

1002.9 

100 

10,0 

53,0 

1001.7 
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7.0 CONCLUSIONS 


Development of an analysis package for the determination of noise genera- 
ted by a moving source with respect to a stationary observer has been accom- 
plished. The procedures outlined in this document when applied to flyover 
data yield a static equivalent noise field with a high degree of statistical 
accuracy. Its utilization for a static/flight comparisons will aid in an 
understanding of forward speed effects on aircraft flyover noise. 



8.Q SYMBOLS 


X position of tho aircraft rnlative to thr microphone 

y position of the aircraft relative to the microphone 

2 position of the aircraft relative to the microphone 

time difference between microphone positions 

distance between microphone i and reference microphone 

average velocity of the aircraft 

directivity angle of acoustic source 

number of Fourier transforms averaged 

number of microphones 

M*N (total numb-'r of averages) 

resulta- power spectral density 

power spectral density of microphone i 

number of degrees of freedom 

error of estimation 


normalized random error 



> * t 




n 

ASPLi 

RE 

ASPL2 

a 

“CL 

“rot 

“vib,N 

“vib,0 

T 
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data window weighting facto*' 

freguency •'^f son»^co 

frequenev at the observer 

Mach number of the source 

pressure of the source 

pressure at the observer 

exponent In convective amplification equation 

difference in sound pressure level due to convective 
amplification 

distance between acoustic source and observer 

difference in sound pressure level due to the inverse square 
law 

absorption coefficient 
classical absorption coefficient 
rotational relaxation coefficient 
vibrational relaxation of nitrogen coefficient 
vibrational relaxation of oxygen coefficient 
temperature 
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> a • 


Kf U '»«5' 


RH 

f 

*^sat 

H 

^r,0 

a(f) 

ASPLg 


G 

k 


Ar 


f 

c 

ASPL4 


SPL 


fkc 


SPL 


fkm 


SPL 


fkB 


relative humidity 
frequency 

pressure of saturation 

absolute humidity in percent 

relaxation frequency of oxygen 

relaxation frequency of nitrogen 

absorption frequency at frequency f 

difference in sound pressure level due to atmospheric 

absorption 

depth of layer i 

ground factor 

wave number 

difference betiveen direct and reflected path 
central frequency 

difference in sound pressure level due to ground impedance 
corrected SPL at frequency f|^ 
measured SPL at frequency 
background SPL at frequency fj^ 
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signal -to-noise ratio at fj^ 

SR sample rate 

BW bandwidth or frequency resolution 

NPTS number of spectral output points 

t|, time increment of assumed local stationarity 

NSHIFT. number of points corresponding to At. 

p2(f) average power at frequency f 

REF reference for dB conversion 

SPL,(f) average sound pressure level at frequency f 

a 
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9.0 APPENDICES 


APPENDIX A. - Program to Shift Microphone Data 


fR0CI»An OUIISH (INPUT, OUTPUT. T«PE1.T*PE2/T*PE3,T(»PE4,T<^S,T*PE6, 
1 TAP€7.TPPEB,T«PEC.T<*P£'ie,TPPEn,T*PE18, 

Z T«PE13,TPPE14-1»#>UT) 

OUIBSH SHIFTS HICBOPHONE DPTrt FROD OUlB DATA. IF N IS THE 
niCROPHOHE NuriBER EACH CHAHHEl GETS SHIFTED FORUARO BV 

TAUmD(H- 1 1/UEDlSP POINTS UHCRE D 15 THE rilCROPHOHC 
SPACINC.UEU IS THE rtlRCS AFT'S SPEED AND SR IS THE SAflPUE RATE 

OF digitization. 

<irOR PUH^Sa, TAU is in INCRCNENTS of 9375.) 

D. -'>R1DLEV 

a/aa/80 

OVIBSH MAINTAINED BY SYSTEM DEUELOPMENT CORPORATION 


DIMENSION DATM<1875.1E),NAMES<ia),IUNlTS(ia),IHDR{B).IST(ie) 
DIMENSION JHDR(H) 

NAMEUST.'INPUT/'IST.NPTSS 

DATA IST/5,lO.I5,aO,H5,30,35.<<O,^S,50/,NPTS5/l875^, 

1 JHDR^<tJ0H / 

IBLNF-IOH 

CONTINUE 


READ SKIP FACTORS FOR CHANNELS 


read (M. INPUT) 

IF (EOF(M>> EC2.203 

202 STOP ’NO NAMELIST FOUND* 

203 CONTINUE 

READ (H.501O) JHOR 
5010 FORMAT(4A10) 

IF (HPTS5 .GT. 1875) STOP 


•NPTS5 CAN NOT BE CT 1875' 


READ FROM TAPEl ALL CHANNELS AND PLACE ON FILES 2 THROUGH 12 

NPTS5 IS J/E OF THE NO. OF POINTS TO SKIP 

1ST IS THE ARRAY OF piUtTlPLES OF NPTS5 FOR EACH CHANNEL 


READ (1) I5N.NCHAN,(NAME5(I),I«l.NCHAN),<IUNITSm,l*l,HCHAN), 
I (IHDRCI ),I*1.8) 

IF (EOFtl )) 5.9 

5 STOP • RECHECK -TAPEl IS INPUT TAPE- SHOULD HAUE 10 DATA CH'S' 

9 IF (HCHAN.CT.12) CO TO 5 
Ll'0 tll*0 

IJ.0 IUK>0 
1 DO 10 K-1.NPT55 

READ (1) (DATA(K.I ).I-1.NCHAH) 

IF (E0F<D) 8,10 

6 LSAU-K-1 
L-LSAU 

IFCL.EO.0) CO TO 21 
Ll-1 

GO TO 11 

10 CONTINUE 
L.NPTSS 

11 MRITE (2) <DATA(K,l >,K-1,L) 

URITE (2) <DATA<K.2).K-l,L) 

IJ-IJ+1 

DO 20 J.3.NCHAN 

IF (lST(J-e>.CE.IJ) GO TO 20 
WRITE U) (DATAfK,J),K-l,L) 

IF (J.EQ.NCNAN) IJK>IJKtl 
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45 

S« 

51 

e« 

IM 


COmiMUE 

IF (Li.Ea.«) 60 TO 1 
DO M J>3,NCHAN 
ElHOmE J 
EMDMLE J 
REMIND J 
CONTINUE 
CRU CVICTtl) 

IF (JHDR(l) .EO. IIIjNK) 00 TO 31 
DO 3B INT»1,4 

IHD«<1HT>*JHDR(IHT) 

CONTINUE 

URITE (13) ISN,NCMRN,(rW«ESa),I-l,NC«RN),aUHITS(l),X-l,NCM(WO, 
1 ;<1HDR(I), 1-1,8) 

URITE FILES 8 TWROUGH IP. OH TW»C1 

IJI-* Wa-NPTSS 

RERD <3><DRTR<K,1),K-1,HPTSS) 

READ (8) <DRTR(K.8),)C-1,HPTSS) 

1J1-IJ141 

IFtlJl.EQ.IJK) KL-LSRU 
DO 5« J-3,HCHRH 

READ (J) (DATR(IC,J),K-1,KL) 

IF (EL.EO.LSAU) L-ICL 
IFOCL.EQ.LSAU) Il-l 
IFtKL.EO.LSRV) CO TO 51 
IF (E0F(J))45,5t 
Il-l 
L-LSAU 
GO TO 51 
CONTINUE 
L-NPTS5 

IF (L.E0.8) CO TO IM 
DO 50 I-l/L 

URITE (13) {DRTft(I,J),J-l,NCNAH) 

CONTINUE 

IFai.EO.8) CO TO 48 
EMDfILE 13 fENDFILE 13 
STOP • SUCCESSFUL RUN • 

END 
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APPENDIX B, - Acoustic Analysis Proqrarn 


The Acoustics Analysis Proqraiti is a tinin sonies analysis program main- 
tained by 


System Development Corporation 
3217 North Armi stead Avenue 
Hampton, Virginia 23656 


It utilizes the Cooley-Tukey algorithm for the Fourier Transform which con- 
verts time domain data to frequency domain data. Power spectral densities are 
calculated by the direct method, i.e., the power spectf'al density is propor- 
tional to the Fourier Transform of the data squared. The sound pressure level 
is simply the power spectral density converted to decibels. 

The program also has the capability to calculate, print, and plot one- 
third octave spectra, atito corral ations, cross spectral densities, cross cor- 
relations, coherence functions, and transfer functions. It can also retain 
most of these functions {or further calculations. It is a lengthy program, as 
is its input parameter list. It also is system dependent and will not be 
presented here. 
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APPENDIX C, - Program for Erisernble Averaging 


PfiOCR«H flWOATfl < I NPUT, OUTPUT, Tape 1 , TAPE?, TAPg 1 6 , TAPES* IHPUT, 

, TAP£6*0UTPwT, TAPES, TAPE4' 

PROOPAA AAHDATA 3S A PAC<nCE FC« HANIPUTATINS DATA PROtl A 
rifT FILE. IT eOMSISTS CF BASICALLY TUO Rf.UTlHESS OME TO STAC*‘ 
riLEi? and OHE to PEFCRit the basic PATHEmATICAL ruHCTJOMS. JTS 
CREATION WAS NECESSARY FOR NANIPULATINQ FILES CREATED BY the 
ACOUSTICS ANALYSIS PROORANt aaP^UN*3'»SDVIiC ) AND THE SICNAL ANALYSIS 
procranineusap uN-s-tSMSo; . but jt may bc usfd on any TIFT nu. 
THE stacking routine, mCuEyER, IS CNE PPRE APPLICABLE TO AAP AND 

hewsap created files. 


DOREEN GRIDLEV 

JANUARY 1580 


NAHPATA NAINTAINE: By SySTF.t CEUELCPPENT CCRPORATJON 



tttx 

NANELIST flNPljTB ZfU 

variable 

DEFAULT 

DESCRIPTION 


ttttmii 

fmtiszxixrxZH 

lOPT 

INTEGER 

0- BO NOT STACK FT-Ef 
1* STACK ALL SERIAL HUTfiERS 
5* STACK SELECTED SERIAL NUNBERS 


» 

JSN 

array 

SERIAL WPBFRD TO BE STACKEDIUSED 


J©0»o 

UHEN JOPT*K) 

EKANOLF.! JSN.1,3,4 InPLIEA 
That SERIAL NuNDERS I. 3, AND 4 
WILL BE stacked 
IF J6N.f,«>5, A constant UIU. BE 
UR’TTEN TO CHANNEL I. 

CONS'' 

REAL 

constant to be ENFLOVED UHEN 


1 

JSN*9S3 

IFUNC 

IN'-EGER 

0-NO ACTION 


0 

l-ADD CHANNELS 

2*5UBTRACT CHANNELS 

3*hUl.TIPLv CHANNELS 

«1*DIUIDE channels 

S"AUERAGE CHANNELS 

6-CHl -SQUARED EQUIVALENCE TES" 

EXAMPLE! 10PT*e.JSN*i ,2.^,IFUNC*i 
SERIALS 1,5, AND 4 UILL BE 
stacked as CHANNELS 1,S,AJ<D 3 
RESPECT I UELV. CHANNEL 3 UILL 
BE subtracted FROn CHANNEL 1 
A-: UILL CHANNEL 5. THEREFORE, 
';FE RESULT, l-e-3 


KSN 

INTEGER 

NEU SERIAL number or FIRST FILE 


0 

REGAINED (SUBSEQUENT FILES. IF ANY, 
WILL BE IN sequential ORDER PROM 
KSN) 

ISPL 

INTEGER 

©•INPUT 16 NOT IN BB 

e 

1*1NPUT IS IN DB 

OSPL 

INTEGER 

O-OUTPUT NOT TO BE IN DB 


0 

1 "OUTPUT TO BE IN DB 
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c ut 

C lit 

c tit 

j 

C tit 

C s« 
C Itt 
C ttt 


lee 


liO 


lee 

c «*t 

139 

140 
1^0 

200 
C sit 


210 


6000 

220 

250 

270 


» ^ 


ORISINAL PAGE tS 
nc onnn OUAUTi 


•i 

I 

r 


REF 

IIMR 


REAL CONVERSION FACTOR FOR JS^L AND OSPL 

o.oooe 

INTEGER FLAG TO CALCULATE ASSVNPTOTIC AND 

0 estimated variances UHEH /iUERAQINS 

0»D0 NO calculate 
1 'CALCULATE 

stsiststttis 


DIMENSION JNAMEStiei l.JUHITSdOl I, JSN( 100 ), DATA( 101 J.NAMESdOl >, 
lUHITSdOl J.HDRtBl 
REAL OSPL*JHDR(S) 

DATA lOf>T,JSN,OON5T,IFUNC,KGHd6PL,OSPL,REF/l,lO0*O»8*O#O,2»O» 
0.0002/, lUAR'O* 

NATtE IST/INPUT/IOPT, J5H, CONST, 1FUMC,K6N,I5PL,05PL, REF, IVAR 

READ NAMELIST IHPU’' 

READ INPUT 
IF {EOPCSn 9999,2 

THIS SECTION IS FOR STACKING SERIAL NUMBERS 

THAT CONTAIN ONE CKANMEL OTHER THAN TIMEfOR FREOUENCV) 


K*1 

URITF (6, INPUT) 

IFdOPT.EO.O) GO TO 1000 
N*1 

REUIHD I 

READ d ) IEN,NCMAN, (NAMEfH I) , I > 1 ,HCHAN ) , < lUHiTSi 1 J, I • 1 .NCHAH ), 
(HDRtI),I*l,B> 

IF CEOrtD) 900,110 

CONTINUE 

IJK'O 

IF dOPT.EQ.l ) GO TO gOO 
KN-H-1 

IF (KN.EO.O) CO TO 150 
DO 120 1*1, KH 

IF (ISH.EO.JSNd )) CO TO 130 
CO TO 150 

IF NEEDED, CHECK TO SEE IF ISN IS A DESIRED SERIAL NUMBER 

READ d) 

IF (EOFCin 100,130 
IF CJSN(N).E0.999» GO tq 200 
IF dSN-JSN<Hn 130,200 
K-K + 1 

K IS the Number of isns fouhd+i 

IF (K.NE.2) GO TO 450 
IF <JSNIN).EQ.999> GO TO 300 
DO 210 I»l,2 

JNAMES(I)*NA«£Sd» 

JUNITSCI )-IUNIT5d) 

ENCODE de,G0&O,JUNITS(2)) ISN 
FORMAT (110) 

DO 220 t«l,8 

JHDRd)'HOR(I> 

READ d) (DATA<n,l«l,2) 

IF (EOFd)) 600,270 
IJK-IJK+1 
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IF nJ<.E0.a> TI«E*DflTrt(l) 

IF (JSM«HJ,NE.P»9) CO TO <<00 
DOTom-0. 
pOTfl(2).COHST 
IT0«“ 

IF {MOP(»C,aJ.EO.O)JTO<<^ 

UPITE (ITO) {DATA(n,I.J,31 
CO TO aso 

IF (JCmH).NE.98»} CO TO SOO 
JMWIESC K ) • tOHCONSTAMT 
JUMITSaO-lOH 999 
SDATrt. CONST 
CO TO 5U 

£MCOr€llO,e001,JUHITS(Kn ISN 
6001 FORMAT! 1 10) 

JHAnE5tK).NAMES<g) 

READ 11) (0ATA<L),U<>1«2) 

IF !EOF!l))600,510 

ISDATA-DATfti?! 1 flPDAfrt«PftTAl 1 ) 

IF ajU.NE.g) CO TO Ell 
IF iDATACn.EO. TIME) CO TO PU 
IF t jeH5N“l).EQ,99«) fit) TO 511 

PRIHTt, 'SERIAL NdMBFR MEN,' TIMES DO NOT AGREE WITH PREgioUS St 
•RIAL* 

STOP 'STACKINC NON-COMPATADLE 6ERIALS' 

J<*K 1 

IF (JSH(N),NE,999 .AND. UK .EQ. B) TinE-DATA'l) 

I TO- 7 

IF (MOOt);,g).EO.0)lTO-8 
JTO'8 

IF (ITO.EQ.8) JTO-7 
READ <JTO) (DATA(L).L-1,JK) 

IF (EOFlJTOn GO0,R20 

)<i4DATA 

IF (JSN(H).NE.S99> DATAf))-FDATA 
URITEtITO) (DATA(L5,L-1.K) 

IF U5NCN).NE.999)G0 TO EOl 
SDATA-CONST 

CO TO sia 

ENDFILE ITO 
ENDFlLE ITO 
IF (K.NE.2) REWIND JTO 
REWIND ITO 
H-N+1 

IF tIOPT.EO.n 00 TO 100 
IF (JSNtNl.HE.O) CO TO 100 
DO 67S I-l.K 

lUNITStn-JUNITStn 
NAMEStn-JMAMEStl) 

DO 680 I-UK 

HDRtn-JHDRtl > 

KSN-KSNtl 
liCHAN-K 
KTO-4 

WRITE (KTO) x:5N,MCHAH,(NAr5EStn»I-l.NCKAN!,(IlJHlTS(n,l-l,NCHAH), 

. (HDR(l), 1-1.8) 

read (ITO) (PATA(I),I-1,NCHAN) 

IF (EOFtlTO)) 710.730 
EHDFILE KTO 

PRINT *. • WAUUES ARE STAC#:ED ON TAPE4* 


300 

^00 


•<60 


600 


601 

510 


5U 

61S 


550 


600 


670 

676 

680 


700 

710 


(TCUIND ITO 
PEUINP JTO 
CO TO 1000 

720 WRITE (KTO) (DflTRd ),1»1,«CHAN) 

CO TO 700 

000 IF IJOPT.EO.n CO TO 670 

STOP • recheck 5ER10LS— one IS NOT FOUND* 

1000 IF tlFUNC.NE.O) CO TO tOOl 
C 

C TO CONTINUE TO STOCK SERIOLS ON TOPE-I 
C 

READ INPUT 
IF (EOFiei) 1002,2 
1002 ENDFILE KTO 

STOP 'REOUESTEO CHANrCLS STACKED OH TAPE^* 

1001 IJ-0 

IF (lUAR.HE.OI URITE <6,5020) 

5020 FORflAT /,10X,IFRE0UENCV *,5X,*AUC FWER »,EX,» UARA », 

1 SX,R UARE*,//) 

REWIND N 

READ t tt 1 1 5N , f t HAfKS ? !) , I • 1 , NCHAN ) , < I UN 1 Tg { | > , I ^ 1 , NON AH ) , 

, ■ «HDR(1),1*1,B» 

IF CEOF(ri)) lose, 1060 
1050 STOP ’TAPE-t IS EMPTV 
1060 NCH-NCNAN 
NCHAN.2 
IJKL»0 

1070 READ (M) (DATAUl.l.l.NCIO 
IF C£OF(M)) 2000, 1080 
1080 IJKL-lJKL+1 

1F<ISPL.E0.1.AHD.IJKE.0T.7) CAU FDOCON'.DA’^A,NCH,REF ) •IJ«IJ+1 
GO TO (1100,1200,1300.1400,1500.1700) IFUNC 
STOP * IFUNC CODE NOT DEFINED PROPERIH* 

C SXt 

C ««» ADDING CHANNELS A+B+C+D+... 

C «»» 

1100 NAnES(2)-10HADDED CHS 
ADATA-0 

DO 1125 1-2. NCH 

ADATA»ADATA+DATA(I) 

112S CONTINUE 
GO TO 1600 

C tX!t 

CSS* SUBTRACTING CHANNELS A-B-C-D-E-F- . . . 

C SSS 

1200 NAnES<2)»10HSUBTR CHS 
ADATA*DATA(2) 

DO 122S I *3, NCH 

ADATA-ADATA-DATACI) 
laas CONTINUE 
Oo TO 1600 

c ss* 

C ss* MJLTIPLVIHG CHANNELS ASBsCSD*... 

C SSS 

1300 NARE5(2)-ief«JLTI CHS 
A DATA *1 

DO 132S I-a.NCH 

ADATA-ADATA»DATA(1) 

132S CONTINUE 


* - #,7V.. V' " •- - 

, .ir’- N V •• ’ 

OP pool? 


fV*’'’ 

!* i ' U 


c *t» 
c <tt 

C St* 
HM 


14c5 

C tt* 

c stt 
c tt« 
i5&e 


15BS 


1550 


5010 

1700 
C lit 
C ttt 

C »H» 


1600 


aeoe 


8999 

2000 

c its 


SOM 

5001 


GO TO 1600 

DiyiDlMO CHAMM6LS 0/B/C/D/. . . 

HAnes^sj-ioHDiuiD chs 

ADAT0*DATA(a> 

DO 1485 1-3,NCH 

ADATA-ADATA/'DflTAm 
CONTINUE 
CO TO 1600 

AUERACINC CHANNELS (A+B*0+D+ . . . )/(NCH-l ) 

HA«ES(aMOHAUGED CHS 
sun-0 

DO ISaS I-2.NCH 

SUM-SUN+DATAU) 

CONTINUE 

ADATA*SUn/(NCH-l ) 

IF UUAR .EQ. 0) GO TO 1600 
MARE *0 * 

MARA-ADATAtira^(NCH-l ) 

DO 1550 I-a,NCH 

MARE-VARE+ C t DATAO )SS2 )-{ADATAt*e ) ) 

CONTINUE 

MAR£-MARE/<(NCH-nt»a) 

URITE (6,5010) DATA(1),ADATA,UARA,UAR£ 
FORHAT(l0X,3(Eia.5,5X),Eia.5) 

GO TO 1600 
COHTirflJE 

CHI-SQUARED EOUIUALENCE TEST 
IF (IJ.EQ.l ) SDATA-0 

SDATA-SDATA+(ALOG10(ABS(DATAf3)/DATA(3 ) )>)*»a 
GO TO 1070 

DATA(2)-ADATA OIFClJ.EO.l) KSN-KSN+1 
IF (OSPl.EQ.l.AND.IJta.CT.?) CALL TDBCON(DATA(a ),REn 
IF (IJ.EQ,1) URITE (10) KSN.NCHAN,(HAMES(I ),l-l,HCHAN), 
(lUNITSCI ),I-1,NCHAM),(HDR(1 ),I-l,8) 

WITS (10) (DATA(K),K-l,a) 

GO TO 1070 

IF (IFUNC.EQ.6) GO TO 3000 

EHDfILE 10 

REuiND M 

CO TO 1 

EHDFILE 10 

STOP 'END OF PROCRAn* 

CONTINUE 

COmiTE CHI -SQUARE 
CHIS<J-IJ»SDATA/‘^. 

DF-FLOAT(IJ) 

PRINTS, "DF- ■ , DF, 'SDATA- • ,SDATA, "CHI GO* " ,CHISQ 
(NILL nDCH(CHISO.Df,P,IER) 

IF (lER.EQ.jaS .OR. IER.E0.34) PRINTS,* ERROR IN USING CHI* 
P-(1-P)li00. 

WITE (6,5000) 

F0RnAT(/'/,15X.» CHI-SOWiR£ ECOIMALENCE TESTS.-//) 

URITE(6,S«01) NAI1ES(a),NA«:S(3),P 

FORnAKSX.S THE PROBABILITY THAT CHANNEL 8,A10,S AND CHANNEL S, 
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.Alt./.GX.t «RE STflTISTICALLV EOUIVAUffT IS «.F10.5,S PERCENT*) 
END 

SUSROUTINC EDBCONC DATA, NCH, REF) 

DIMENSION DATACiei) 

DO 1 I-a,NCH 

DATAd >•(!«, **<DflTA<l)/ie. ))*REF»REF 
1 CONTINUE 
RETURN 
END 

SUBROUTINE TDBCON<DAT,R£F ) 

IF (DAT) a, I 

a DAT-I0.*ALOCl6(ABSfDAT)/R£F*REF) 

1 CONTINUE 
RETURN 
END 


* 
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APPENDIX D 


ORIGINAI- P&GE ]3 
OF POOR QUAus"! 


. - Program for Correcting Sound Pressure Levels 


PROfiBAn PLVOUER I INPUT, OUTPUT, TAPES* INPUT, TAPE6*0UTPUT, 

1 TAPEl,TAPEa,TAPE10) 

PROCRAfI FLVOUER IS TO CORRECT OUIB FLVOUER SOUND PRESSURE 
LEUEL DATA GENERATED BV THE ACOUSTICS ANALYSIS PROCRA«( AAP ) . 
THE propagation path corrections ARE! 

1> INSTRUNENTATION corrections— PRESSURE RESPONSE, 
diffraction, and windscreen corrections 

2) conoectiue anplification 

3) INUERSE SQUARE LAU 

A) ATHOSPHERIC absorption 
S' ground irpedance 

S) DOPPLER FREOUENCV SHIFT 
FILES! 

1 ) INPUT FILE IS TAPEI 

2! OUTPUT FILE FOR CORRECTED SPL'S IS TAPEI0 

FLVOUER NAINTAINED BV SVSTF.H DEI/ELOPflENT CORPORATION 

the FIRST INPUT EXPECTED IS A SIXI6) DIGIT CODE BEGINNING 
IN COLUNM ONE(l). THE ORDER OF THE CODE IS THE ORDER OF 
THE ABOUE CORRECTIONS. 

0 IMPLIES DO NOT CORRECT 

1 IMPLIES CORRECT 

EX.li 000000 NO CORRECTIONS 

EX. 2! 101000 CORRECT FOR INSTRUnEHTATlON AND INU. SO. 
A NAMELIST IINPUT# IS EXPECTED NEX’. 

NAMELIST SINPUTS 


UARIABLE t TVFE AND x DESCRIPTIOf 
I DEFAULT X 

«*ItSXXX»XXXXXXXXIXX'X 


ICODE INTEGER FLAG FOR WHICH THETAS ON 7APE1 
0 TO CONSIDER 

0*A1.L thetas ON TAPEI 
1 -SELECTED THETAS ACCORDING 
TO ARRAV THTp 


THTA REAL ARRAV ARRAV OF THETAS TO BE CONSIDERED WHEN 
10X0,0 ICODE-1, 

NOTE! THETA ON FILE TAPEI UJLL BE 
ACOfiPTED IF IT IS UITHIN ♦.! 

OR -.1 Of THTA 


N 


INTEGER FLAG FOR TYPE OF ACOUSTIC SOURCE 

e DISTRJDUTIONtfOR CONU, AMP.) 

e-nONOPOLE 
I-DXPOLE 
E»QUADRAPOLE 


STRCR 


REAL ARRAV 

sees -999, 


ARRAV OF UALUES 
CORRECTIONS 

STRCR a, n ARE 
STRCR (1.21 ARE 
STRCR (1,3) ARE 
STRCR (I. A) are 


FOR INSTRUMENTATION 

THE FREQUENCIES 
DTSPL FOR P. RESPONSE 
DTSPL FOR diffraction 
DTSPL FOR WINDSCREEN 


1 
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or* pcb»: 


HGT 

REAL 

30. 

DISTANCE FROM GROUND TO fllCROPMOfiEtFEET ) 

REP 

REAL 

15. 

REFERENCE DISTANCE FOR INVERSE SOUARE 
LAU CORRECTION 

P 

REAL array 
35T-999. 

MEASURED ATMOSPHERIC PRESSUE <NBARS) 
CORRESPONDING TO ARRAY ALT 

T 

REAL ARRAY 
2RI-999. 

NEASURCD ATMOSPHERIC TEMPERATURE (DEG C> 
CORRESPOND INC. TO ARRAY ALT 
NOTE! THE TEMPERATURE ARRAY MUST 

BE INPUT WITH EACH SUCCESSIVE 
NAMELX.S7 IF APPLYING CORREC- 
TIONS 2, -1,5, OR 6 

ALT 

REAL ARRAY 
25*-999. 

MEASURED ALTITUDE (METERS) 

input is to be in DECREASING ORDER 
AND IN EQUAL INCREMENTS 

RH 

REAL ARRAY 
251-999. 

MEASURED RELATIUE HUMIDITY (PERCENT) 
CORRESPONDING TO ARRAY ALT 

lUAR 

INTEGER 

0 

FLAG TO PRINT OUT RESULTS OF SPL 
AFTER F.ACm correction IS HADE 
0"D0 NOT PRINT 
1 "PRINT 

FJ 

REAL 

0.0 

INITIAL frequency FOR CORRECTIONS 
TO BE APPLIED 

ra 

REAL 

50000. 

FINAL FREQUENCY FOR CORRECTIONS TO 
BE APPLIED 

SR 

REAL 

50000. 

SAMPLE RATE OF DIGITIZATION 

NPTS 

INTEGER 

SIS 

NUMBER OF POINTS PER BLOCK 
(UARlABI.E NREAD IN AAP) 

NBLKS 

INTEGER 

5 

NUMBER OF BLOCKS FOR THE SPL USED 
IN AAP 

NMICS 

INTEGER 

8 

NUMBER OF MICROPHONES AVERAGED 
IN MAHDATB AFTER AAP 

ITAPE 

INTEGER 

0 

FLAG FOR WRITING (AND PRINTING) ONLY 
THOSE FREQUENCIES BETWEEN FI AND F2 

0- URITE ALL FREQUENCIES 

1- URlTE FREQUENCIES BETUEEN FI AND Ffi 

IPRINT 

INTEGER 

FLAG FOR PRINTING WHEN IVAR»i 


©"PRINT DELTA SPL UALUES 
l"PRUr CORRESPONDING SPL WLUES 


FOR rt. mJELI.ER 
Ik-81 
D.Cr<rDLEY 


55 




C (THIS IS REUISED WERSIOM Of ft PROCRRM DEUELOPtO 5-8®) 

C 

DIMENSIOH THTA(10),6TRCn(2Oe,4>,P(gS>,T<25),flLTC25),RH(2S),IFHC(6> 
1 ,MTfltE5,NAriE5(2),IUNITSl8),IHDR<8) 

HiVieUST /INPUT/ ICODE,T^^Tn,H,ST^CR.HGT,r?£F,P,T,PLT.RH, 

1 Fl,Fa,IUftR,lTftPE,IPRIHT^SR,NPTS,NBUS,Mf1ICS 

WTR ICODE,THTfl,N,STRCR/0, 10*0.0.0, 8005-999. /*HGT/30./, 

1 P,T,PLT,RH/100»-999./,IURR. FI, F2/0.0., 50000./, 

a NBLKS,NMICS,ITW>E/S.8,O/,IPRrNT/O/,SR,NPTS/50000.,5ia/ 

MTfl CORR, DTSPL, D2SPL, D3SPI., D-1SPL , DSSPL/1 .,5*0.0/, REF/ IS . / 

C 

C READ PRRflV IFNC FROfl INPUT 
C 

1 READ (5,5111) (1FNC(I),X*1,6) 

5111 F0RMAT(6I1 ) 

IF CEOFcsn 1000, a 

C INITIALIZE UARIABLES 

C 

a PI'3.H159265 
PO-1013.35 
T0»a93.1S 
A-0,01 
C 

C READ NAMELIST *INPUT S 

C 

3 READ (5, INPUT) 

IF (E0FC5J) 1000, -1 

4 CONTINUE 
C 

0 URITE INPUT INFORMATION TO OUTPUT FILE 

C 

URITE (6,5322) 

52aa FORMAT! 1H1,//,5X,*1FNC • CODE FOR CORRECTIONS!© MEANS NOT APPLIED* 
1,* 1 MEANS APPLIED)*/, SX,*IN THIS ORDER:*, /, lOX, 

2 INSTRUMENTATION CORRECTIONS*, /,10X, 

3 ICONUECTIUE AMPLIFICATION*,/, lOX, 

4 5INUERSE SQUARE LAW*,/, lOX.XATMOSPHERIC ABSORPTION*,/, 10X, 

5 »GROUND IMPEDANCE*,/, 10X,*D0PPLER FREOUENCV SHIFT*,/) 

URITE (6,5333) ( IFNC! 1 ), 1» 1 ,6 ) 

5333 FORMAT! 10X,»IFNC" 8,GU) 

URITE (6, INPUT) 

C 

C READ HEADER RECORD. AND TR, RE, THETA, X.UEL.Z, AND V (THE FIRST 7 

C data records 70 DETERMINE IF THE FILE IS ONE REQUESTED BV ICODE 

C AND THTA 
C 

KCNT-O 

le READ (1) ISN,NC)IAN,!NA«E 5!1 >,I*l.NCHAN),aUNITS(n,I'l,NCHAN), 

1 (IHDRd ). 1*1,8) 

IF (EOF!!)) 900,20 
e« IJ*0 

as READ (1) (DATA!! ),I*1,NCHAM) 

IF (EOF(l)) 30,40 

30 STOP 'END OF FILE BEFORE RADAR INFORMATION IS READ’ 

40 IJ»IJ+1 

IF (IJ .EO. 1) TR-DATAO) 

IF (IJ ,EQ. a) R£*DATA(3) 

IF (IJ .EQ. 3)T)«TA-DATA(a) 

IF (IJ ,EQ. 4) XRAD-DATA(a) 
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IP (IJ .EO. B) yEL»DATfl(H) 

IF (IJ .EQ, 6 ) Z£«MTn(S) 

IF <1J .EO. 7) VRAD-MTrt(2) 

IF (IJ .NE. 7) M TO as 

too IF (ICODE .EO. 0) CO TO ISO 
DO no 1-1. 10 

IF (THETA .CE.(THTA<n-0.1). AND. THETA. LE.(THTAm+0.n) 

1 CO TO ISO 
no COMT INLIE 

OTHERWISE CALL SKIPFF 

CALL SKIPFF(5LTAPE1,1 > 

CO TO 10 
ISO KCNT-KCHT+1 

WRITE HEADER RECORD AMD RADAR iMFORtlATIOM TO TAPE 10 

WRITE (10) ISM,NCHAM,(MAfieS(I >,I*l,MCHAN),<IUNnS(IM«l,NCHAN), 
1 (IHDRII), 1-1,8) 

DUM--993. 

WRITE (10) DUM,TR 
DUM*-990. 

WRITE (10) DUH.RE 
bun ■•'■997. 

WRITE (10) DUn, THETA 

Dun--i:!9G. 

WRITE (10) DUn.XRAD 
DUn--9S>5. 

WRITE (10) DUn.OEL 
DUM--99.1. 

WRITE (10) DUn.ZE 
DUn--993. 

WRITE (10) DUn.YRAD 
IF (KCNT .CT. 1 ) GO TO 180 

ADD 373.15 TO TEOP ARRAY 

DO 135 I-l.SS 

135 IF (T(I) .ME. -999, )TtI)'T(I)+273.l5 

COHPUTE AVERAGE TEMPERATURE FOR AVERAGE C FOR DF,ANb CA 

AT-0. 

ICNT-0 

I-O 

140 i-m 

IF (T(l) .EQ. -999. .OR. 1 .CT. SS ) C.0 TO 150 
1F((ALT(I) .LT. (HGT*,3048)) .OR. (ALT(I) .CT, (2EI.304B))) 

1 CO TO 140 
AT-AT + T(I> 

ICNT-ICNT+1 
CO TO 140 

150 AT-AT/(FLOAT(ICMT)) 

TA-AT/TO 

CAVO-(343.S3*SORT(TA>)/6.304B 
COnPUTE DELTA HEIGHT IN ALTITUDE 
DHGt-ALT(l>-ALTCa) 
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COftPUTE C CLOSEST TO CROOHO FOR GROUND IftPEDONCE 

0L70°SD999. 

DO JfiO I-UE5 

7F (ALTCI) ,EQ. "999.) 60 TO 170 
IF <ALT(I) ,CE. OLTO) GO TO d60 
0LTn«ALTU ) 

ID-r 

ISO CONTINUE 

170 CGl*(343.a3«S0RT(T(ID)/T0) )/O.30')8 

CONFUTE UftRIODLES FOR CORRECTIONS 

CORK IS FOR THE DOPPLER FPEOUENCV SHIFT 
ANCH IS THE PUERftCE riACH HUNBER 
DTSPU IS FOR THE INUERSE SOUARE LfiU 
D2SPL IS FOR CONUSCTIUE ftnPUFICATIOH 
DK IS A NUNESR BASED ON nANWJIDTH 
DR IS THE PATH DIFFERENCE DETUEF.N REFLECTED AND 
DIRECT path 

DKDR IS THE PRODUCT OF DK AND DR 

180 ANCH-UEL/CAUC 

CORR- ( 1 . -ANCH»COS(THETA*PI/180 . ) ) 

DTSPL*30.*ALOG1O<REF/RE ) 

IF (IFNC(3) .EO. 0) D'etre. 0.0 
D2SPU»e0 . S (28N+2 )*AL06‘.0( CORR ) 

IF (IFNC(2J .EO. 0) D25PL-0.0 
DK*(Pl»SR)/<FLOAT<HPTS)»Can 
DR-2.»HOr»SINnHETA*Pl.'180. ) 

OLD NETHOD OF CALUCULATING DR--IHUALID DUE TO CHOICE OF ANCLE 

DR*S0RT(REtRE + <1.»ZE»HCT)“RE 
DKDR'DKJDR 

IF (DKDR .EO. 0) Ul*l. 

IF (DKDR .HE. 0) U1-S1N(DKDR l/DKOR 

URITE TABLE HEADINGS OH OUTPUT IF IUftR*l 

IF (lUAR .NE. 1 > GO TO 139 
WRITE (6,5888) ( IHDRO- ), L“ 1 ,8 i 
5838 FORRAT (lHl,/'//.BAlO) 

WRITE (6,5-)4<n TR,RE.T)-(ETA.KRAD.UEL,ZE,YRAD 
1 ,CAUG, CGI, ARCH, NNICS.NBLKS 

544<1 F0RnAT<//.5)<,8T-»,E12.5,<X,nn->l,E12.5,AN,lTHETA«»,El2,5,<X, 
15)<*S,Eia.S,'5!<,lU.a,EI2.S.4S,sZ«S.Eia.5,'tX,»V-»,E12,5,//, 
A30X.8CAUG- S.ElB.S.^X.SOr.I- X.Eia.B.-JX.SN" «,E12.5.///.2eX, 
aSAOERACE OUER mcRO.^HONES, AND a,!-).* BLOCKS!. //, lOX, 

38C0DE FOR CORRECTIONS RAPEJa,/, ISX.alC-INSTRUNENTATIONt . 15X, 
^aCA-CONUECTIUE ARPLIFICATlPMa,/. I5X.al6»INUERS£ SQUARE LAUI./.ISX, 
E*RA«ATN0SPHERIC ADS0,RPTI(>Na./,15X,aGI. GROUND INPEDAMCE*. 15X, 
GIDF-DOPPLER FREQUEHCV SHIFT!) 

IF (IPRINT .NE. 0) WRITE t6,56SG) 

6GG6 F0RrWT(/'.',4X,a FREOUSHC‘'’a . 

1S(AX,*AUQ SPL,DBa).<X,iJFREGUENCVS,/,3SX,l(IC)»,7X,t(IC,CA)a,7X, 
a»<IC.CA,IS)».^X,a(IC.CA,I6.AA)«.lX,a( all ) 8,5X,t(DF)*,^/) 

IF (IPRINT .EG. 0) WRITE (6.G777) 

5777 FORHAT(.'/.<X.« FRE0UE((CVa,4X, 8W.»G SPL.DB*. 

1 S(4X,8DELTA SPL * ),4X,8FR£GU£NCya4XSAUQ SPL,DBX,/',3SX,«( IC )*,9X, 
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3 l#X,»FIMAU,//> 

13C COtfTINUa 

READ AND ApPLV REOUCSTED CORRECTIONS 

aOO READ (1) F,SPL 

IF fEOFU)) 500,310 

atO IF tF .G£. FI .AND.F .LE. FB) 00 TO e»5 
IF aTAPE .NC. 0) GO TO 300 
FFF*F 

SPLl*SPLe-SPl.3-SPL4-SPL5-SPL 
IF (IlMR.NE. 1) CO TO 460 
00 TO 455 

ai5 IF (IFMCU) .EO. 0) 00 TO 300 

INSTRUnCNTATION CORRECTIONS 

FF*DDl*DBa*DB3«0.O 

IF CF.GT. STRCRf 1, 1))G0 TO 330 
IND*0 
00 TO 370 
330 DO 360 1* 1,300 

IF (STRCR«I,l) ,EQ. -999.) GO TO 380 
IF (F.LT. STRCRd,!)) CO TO 370 
IND-1 

FF-STRCR(I,1 ) 

DBl-STRCRn,3) 

DB3-STRCR(1,3) 

DB3-STRCR(I,4) 

360 CONTINUE 
370 IND1-IND»1 

RATI0-(F-FF)^(STRCR<IND1,1 )-FF) 

DBl* DBl + RATI0»(STRCR<IHDI,3)-DB1) * DBl*INT(DBltl0.+.5)/10. 
DBS- DBS ♦ RAT10»(STRCR(IND1,3)-DB3J • DB3-INT(DB3*10.+ ,S )/'10 . 
DD3- DB3 + RAT10»(5TRCR<1ND1,4)-DB3) ® DB3-INT<DB3*10.’*. 5)^10. 
380 D5SPL-D61 4D63 4DB3 
SPL-SPL+DSSPL 

300 IF <IFNC(S) .EO, 0) CO TO 310 
CONVECTIUE AMPLIFICATION 
SPL-SPL+DaSPL 

310 IF (IFNCO) .EO. 0) GO TO 330 
IMUERSE SQUARE LAU 
SPL-SPL-DTSPL 

330 IF (IFNC(4) .EQ. 0) GO TO 400 

ATMOSPHERIC ABSORPTION 

D4SPL-e.e 
DO 350 1-1,35 

IF (ALT<I> .EQ. -999.) GO TO 350 

IF ((ALTO) ,LT. (HOT*. 3048)) .OR. (ALT(I) .GT. <2E».3048)>) 

1 GO TO 350 

C-343.33XSQRT(TCI )/TO)BlO0. .'(3.54*13. ) 
SATR»10.7V586*(l.-(S73.ie/Ta ))) -5.O38O3»ALOO10(T( I l/aTS. 16 ) 

1 ♦!, 60474 * 0 . 0001«( 1 .-10.88(“Q.H96938( (T( I )/373 .16 )-l . ) ) ) 
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CF 


rx rr 






a ♦«H2873««,Wtl t<le.t»{^.7695S«<t.-(a73.lS/T<P))l-l.) 

3 -3.21B5BB3 

PR»ie,l»SftT« 

Mun*«H(niPR»P(n/po 

FRl.O«(P<n/PO)*t Z4. +4^ie0»HUn«((0.05+HUfl)/<«.Ml4HUNn) 
FRLH*(Pm/PO)»SORT(TO/T(n>»< S. + 3B0.tHUH* 

1 EXP(>6.14a*(»T<I I 

PUPHA* B.6«e*SO«T(T<I)-'TO)l(F«F*PO/'P<n)» 

1 a.84fio.«»(-ii. ) ♦ a.i9*a.000j»(to/T(n)t(p<n/po)« 

a (aaast./Tdmte. * exp<-bs3b./t(J)) / 

3 <FRLO f(F»F/FRLO>J ♦ B. J.6»O.0OOl*(TO/TC O )» 

4 <P<n/PO)»<335g./Ta))?(33S2./T(n>E 

5 EXP(-335e./T(n>/{FRt.N'» (FSF/FRLN) ) ) 

D4SPL-D4SPL ♦ ALPHAsDHCT/tSlN(THETP*PI/180, )> 

3Sa COHTIHUE 

SPL^SPL+D4SPL 
360 CONTI HUE 

400 IF UFMCCS) .EO. 0) CO TO 410 
C 

C CROONP IMPEDANCE 

C 

ourt»<-cAEa.*Pi»rtDR/CGnK*a. ) 

IF <WJfl at. -STS. 84) H3-0. 

IF (DOtl .GE. -675.84) H2-EXPtDUn) 

D3SPI*10.*ALOC10 <2. ♦e.»U2»C0S<2.*PUF»DR/CGI ).IU) ) 

4CJ3] • 1 

410 IF (IFHC<6) .EQ, 0) GO TO 420 

0 

C DOPPLER FREOUENCV SHIFT 

C 

F-FECORR 
420 CONTINUE 
C 

0 SET UP AND PRINT UALUES IF IUAR*I 
450 IF (lUAR ,N€. 1) GO TO 460 
IF (IFNC(6) .EQ. 0) CORR-I.0 
FFF'F/CORR 

6PL1*SPL+D3SPL-D4SPL+DTSPL-D2SPL-D5SPL 
SPL2*SPL+D3SPL-D4SPL+rT6PL-D2SPL 
SPU^•SPL♦D^SPL-D 4 SPL^DTSPL 
SPL4-SPL*D3SPL-D4SPL 
SPL5*SPL+D3SPL 
465 IF (IPRINT .NE. 0) 

lURITE (6,5555) FFF,SPU,SPLa.SPu3,SPL4,SPU5,SPL,F 
IF (IPPIHT .NC, 0) GO TO 460 
D7SPL— DTSPL 
D6SPL*-D3SPL 

UR1TE(6,SS5S) FFF,SPU,05SPL,D2SPL,D7SPL,D4SPL,08SPL,F,SPL 
SS5E FOPriAT(2X,F10.4,8(4X.FtC.4)) 

46» WRITE (10) F,6PL 
GO TO 200 
500 ENDFILE 10 
CO TO to 
900 REWIND 1 
CO TO 1 

1000 ENDFILE 10 

STOP *END OF PROGRAM* 

END 
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APPENDIX E. - Program for Acoustic nirnctiyities 


PROCy?(Vt SPLTMTC < INPUT, OUTPUT .TPPE5. IW»UT, TAPEfi. OUTPUT, 
1 TAPEJ,TfiPE2, TAPES, TAPEIO) 


PBOCPPn 5PLTHTC IS P PROORArt UHICH GENERPTEO THE 
PADJATIOH PATTERN iSPI, UERtiUS THETA) AT A DESIRED 

rREOuEfi:'! POP AU rites on tapei tvicH are generated 

,6V THE ACOUSTICS ANAtVSlS PROGRAHiAAP ) UJTM PCD FILES 
OVAAP3D,OV'AAP12, AND OUAAPGD APPLIED. 

FOP OUID FORWARD CPEED EFFECTS 
(A, nuELLER) 

D. ORIDLEv (5981 1 

SPLTHTC rtAIHTAlNEP PV SVATEfi PEUEtOPNENT COPPORATICN 
I. FILES 

A. TAPEl IS the PPL Ftl.E CENCRATEp flv aaF.-THE 
FIRST T DATA RECORDS ARE RADAR JHFORNATJON 

B. TAPES IS the radiation PATTERN FILE FOR THE 
BACKGROUND RUN (OMLV USED AS INPUT UHEN 
1D5UE-1 1 

C. TAPES IS the new RADlATimi PATTERN FILE 

D. THPEIO is the new RADIATION PATTF.RN NIHUS THE 
BACKCROUND (TAPES) lONLV OUTPUT UHEN IBSUB*!) 

n. naneust inane* 

A. FREO 15 THE DESIRET FREGUENCV AT UHICH 70 CENEPATE 
THE RADIATION PATTFRT 

B. ERR IS THE bandwidth 0'>m UHICH TO FIND SPL UALUE 
(Range of freo-err «’rf,a*erp is used* 

C. IBSUB IS THE FLAG FOR SUBTRACTING BACKCLROUNC 
(0 • DO HOT subtract 1. SUBTRACT) 

D. 15UHTVP IS THE FLAG FOR CONPUTATIONS 

Ce*Hl CHEST UALUE IN RANGE OHLV.I.SUH OF TUO 
HIGHEST SPL UALUES IN THE RANGE) 


DirtENSION NAnES(2J,lUNlTSc2).1HDP(B),0ATArE),JNAnES(a>,JUNITS(a), 
J THETAB( 50 ),DI«B(EO) 

NftNELlST ,<NAnE/FR£ 0 , ERR , IBSUB . ISUHTVP 

DATA ERRv'aeO/',FREO-'aeO 0 /.IBSUB/O/',ISUrrTVP.' 0 /', 

1 JNAMES/ 10 H theta ,J 0 H SPL 

a JUNITS/lOri DECREES , 50 H DB / 

1 READ ( 5 , NAME) 

IF (E 0 F( 5 >) 1000.5 
5 WRITE ( 6 , NAME) 

K «0 

to READ (S) ISH.HCHAN. (NAMESd ),I-l,NCHAN),<JUNJ 7 S(I ),l*l,NCHAN). 

1 UHDRlI ),I*l, 8 ) 

IF <E 0 F( 1 )) 500 , .ao 
ao K-K+I 

IF (K .EG, 1 ) WRITE ( 9 ) ISN.NCHAN. UNAMES ( I ) , I « 1 .NCHAK), 

I ( JUNITS ( I ) , I • i , NCHAN ) , ( IHDR t I) , I « 1 , Q ) 

J*0 

DDMAX»- 99999 , 

DBinAK«e,o 

30 READ Cn <DATA( 1 ), I»l, NCHAN) 

IF tEOFd )) 35,40 

3 S STOP 'FREOUENCV REQUESTED IS TOO LARGE* 

40 J-J41 
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IF (J .CT. 7) 60 TO SO 
IF CJ ,io. 3) THCTA«DrtTrt(a) 

CO TO 30 

%% IF (DATAd) ,GT.(FRE0 ♦FRR>) 60 TO 6* 

IF IDATAd) ,LT.(FRE0 >ERR>) 00 TO 39 

FIND fWXinun SPL OALUT IN RAHC€ 

IF (ISUNTVP ,E0# B) 60 TO 55 
IF tDATAO) ,LT. DBtfWX) CO TO 30 
IF iDATAO) .LT. DBNAXJ CO TO 52 
DIIHAX-DBAAX li)BNAX*DATA(2> 

GO TO 30 

52 DBINAX*DATA(2) 

CO TO 30 

55 IF {DATA(2) ,CT. D8NAX) DBNAX-DATA(2) 

CO TO 30 

60 IF (ISUfITyP ,E0. 0) GO TO 65 
DBNAX*UO.tt(DBriAX/20. ) 110.0002 
tSBlf1AX-aa.»,l(DBl«AX-^20. ) )»0.OO02 

OBnAX«lS,;ftLOCJ0{<p.SnAXiDBnAH’i'DSifiHXipBifiAX).'< fOoSgl.OOOa)) 

65 IJRITE .O) THETA, DimAX 
READ 1 1 ) 

IF (EOFin) 10,70 
70 CALL S'nPFF(5LTAP.Gl,n 
CO TO 10 
50© EH'DfILE 9 
REUIND i 

IF (IBSUB .EO, 0) GO TO 1 
ENDFILE 9 
REWIND 9 
REUIND 2 

PLACE BKC INTO ARRAY THETA8 AND DBB 

510 RF.AD (21 ISN 

IF CE0F(2>) 520,530 
520 STOP "NO RADAR FILE INFO.* 

530 IJ>0 
540 IJ>IJil 

READ (2) THETABO.n.DBBUJ) 

IF (E0F(2)) 550,540 
550 »CNT*IJ-1 

READ SPL US THETA AND SUBTRACT BXC— ONLY THOSE THETAS 
UlTH A CORRESPONDING TKETAB UILL BE WRITTEN TO TAPEIO. 

60« READ (0) ISN.HCHAM,<HAn£S<l),I*l,NCHAN>,<IUNITC(n.i«l.iKHAN), 

1 (IHDRd ),I-1,8> 

IF (E0F(9)) 1000,610 

610 WRITE (10) ISN.NCHAH.(KAKES(I>,I«l,NCHAN),(lUNITSa),I-l,NCHAN>, 
1 <1HDR(I),I‘1,0) 

620 READ (9) THETA, DB 
IF (£0F(9)) 700,630 
630 L*0 
C40 L>L41 

IF (L .CT. KNT> CO TO 620 

IF (THETAKL) .CT, <THET<VH.) .OR. TKETAB(L) .LT. (TMETA-l.)) 

1 CO TO 640 
650 CONTINUE 
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ir(«DJ-DII(L)) *GT. 3. > CO TO (MO 
IWITE <6,5666 > ISH,THCT6 

666# FORHnT(/,» FOB SERIAL *,I5,I S/W IS < OR • 3 FOR AMCL6 »,riO.<J 
IF (DI>DBI(L) .LE. 0> 00 TO 676 
GO TO 66S 

666 IF ((DB-DBB(L)) .C£. 10.) GO TO 676 

66S DB>OB«10.tAL0010a.-tQ.lt((DBS(U'*DB).^16. )) 

676 URITE (16) nCTA.OB 
CO TO 626 
760 EMDFILE 16 
EHOFILE 18 
CO TO I 

1600 STOP *EMD OF PROGRAM* 

END 
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